Universal scaling of Coulomb drag in graphene layers 



Irene D'Amico 1 '^ and S.G. Sharapov 1 ^] 

1 Institute for Scientific Interchange, via Settimio Severn 65, 1-10133 Torino, Italy 
2 Istituto Nazionale per la Fisica della Materia (INFM), Italy 
(Dated: February 2, 2008) 

We study the Coulomb drag transresistivity between graphene layers employing the finite temperature density 
response function. We analyze the Coulomb coupling between the two layers and show that a universal scaling 
behavior, independent of the interlayer distance can be deduced. We argue that this universal behavior may be 
experimentally observable in a system of two carbon nanotubes of large radius. 
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Coulomb interactions are responsible for a rich variety of by 
phenomena in low-dimensional systems. In particular, there 
are many indications that in semi-metals, such as graphite, the pi\ (T) 

Coulomb forces remain long ranged due to the lack of conven- 
tional screening. This property of graphite has intensified ef- 
forts from both experimental 1 1] and theoretical |2i|3j,|4|] sides 
to understand its electronic properties and interaction-driven 
transitions. However, a direct linear-response transport mea- 
surements of the effect of electron-electron interactions in sin- 
gle isolated sample are usually difficult because for a perfectly 
pure, translationally invariant system, the Coulomb interac- 
tion cannot give rise to resistance. It has been demonstrated 
0] that when two independent electron gases with separate 
electrical contacts are placed in close proximity and current 
is driven through one, the interlayer electron-electron interac- 
tion creates a frictional force which drags a current through 
the other. This so-called Coulomb drag effect 0] is unique in 
that it provides an opportunity to measure directly electron- 
electron interactions. 

The purpose of the present paper is to demonstrate that the 
semi-metal band structure of graphene is responsible for an 
unusual behavior of the Coulomb drag transresistivity. In ad- 
dition our calculations show a universal scaling behavior of 
the Coulomb drag transresistivity, which may be experimen- 
tally observable. 

It is convenient to characterize the Coulomb drag effect 
in terms of the transresistivity p 2 i because in a drag-rate 
measurements p<i\ is directly related to the rate of momen- 
tum transfer from particles in layer 1 to layer 2. How- 
ever, when Kubo formalism is employed ||7|] one arrives at 
the transconductivity a^y- These two quantities are defined 
as P21 = -E^/ii with j% = and ui\ = J2/E1 with 
E2 = 0, where Ei and ji are, respectively, the electric field 
and the current density in layer i. Since the transconduc- 
tivity is caused by a screened interaction between spatially 
separated layers, 0-21 <C <J\\ and p2\ is related to a%\ via 
P21 ~ -o- 2 i/ (o"iiO" 22 ). 

The expression for the dc Coulomb drag transresistivity in 
two dimensional systems, in the RPA approximation, is given E. 



if- 



87r 2 e 2 fc j BTniri2 J 



dKK 



x / du\Vvi{u,K)\ 



2 lmxi(u,K)lmx2(oJ,K) 



(1) 



sinh 



2k B T 



Our purpose is to describe Coulomb drag between two 
graphene layers, i.e. between two planar sheets of carbon 
atoms, so that Eq. Q is purely two dimensional 01 • 

In Eq. (Tfjl Vu(gj,K) = exp(-dK)v(K)/e(uj, K) is the 
screened interlayer interaction, v(K) = 2-Ke 2 / (Ke) is the 
Fourier transform of the Coulomb interaction with the ap- 
propriate low-frequency dielectric constant e, Xii^i K) is the 
noninteracting density-density response function (the "Lind- 
hard function"), m is the total carrier density, d is the inter- 
layer distance, e is the electron charge and i = 1, 2 is the layer 
index. The RPA dielectric function e(u>, K) can be written as 



e{u,K) = [1 - v(K) X i(n,K)][l - v{K)^(Sl,K)] 
- exp(-2dK))v{K) X i(n, K)v{K) X ^, K). 
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We stress that since in contrast to the total current j(K = 
0) = ii(K = 0) + .72 (K = 0), the current in a given layer 
ji(K = 0) is not conserved, the transresistivity Eq. is 
given by fluctuation diagrams (see Refs. [7]), which are simi- 
lar but not identical to the Aslamazov-Larkin diagram known 
from superconductivity |9], and not by a simple bubble as for 
the conventional electrical conductivity. Thus the transresis- 
tivity appears to be sensitive to all kinds of virtual excitations 
that exist in the system whenever lmx(co, K) is nonzero. 

Let us first of all derive the Lindhard function for graphene 
sheets including finite temperature effects. For a single 
graphene layer the Lindhard function may be written as 



x(tfi,K)=x+(ifi,K)+x-(*fi,K) 
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that, setting A = 0, we obtain 
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where n F is the Fermi function, A± 



1 ± 
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k± = k ± K/2, and E± = ^Jh 2 v 2 F k 2 ± + A 2 . Here v F is 
the Fermi velocity, fj, is the chemical potential and A is the 
excitonic gap that may open due to the interaction between 
quasiparticles H. Eqs. 0, 01 and are obtained using 
an effective long-wavelength description of the semimetallic 
energy band structure of a single graphene layer with the dis- 
persion law E(k) 
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simple description in terms of the Dirac equation in 2D (see 
Refs. J3, 0]). In the absence of interactions, in undoped 
graphite, the conduction and valence bands are touching each 
other in the two inequivalent K points. This corresponds to 
choosing \i = A = 0, i.e. in undoped graphite, at T = 
the negative band E = —Hv F k is filled and the positive one 
E = hv F k is empty. Notice that recent measurements of 
quantum magnetic oscillations in graphite 1 10] confirm the 
presence of holes with 2D Dirac-like spectrum E = ±hv F k, 
and this seems to be responsible for the strongly-correlated 
electronic phenomena in this material. The effects of doping 
and opening of the excitonic gap in graphene can be taken 
into account by considering nonzero chemical potential and 
nonzero excitonic gap A, respectively. 

Eqs. 0, and underline the most important difference 
between graphene and the usual two dimensional electron gas 
(see e.g. Refs. IsI fTlll '): while in the latter only the x+ term 
is present (intra-band excitations), in the graphene response 
also the \- term contributes: this reflects the possibility of 
particle-hole excitations with the particle and hole belonging 
to different (positive and negative) branches of the spectrum 
(inter-band excitations). 

Actually when fi = k B T = only the "unusual" \- term 
survives, and one can obtain an exact expression for %(if2 n — > 
uj_+ iO, K) |[l2ll which for A = takes a rather simple form 
1 
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For finite T the response function acquires an additional 
contribution from x+- We find an analytic expression for \+ 
in the limit of small u and K hut finite T and n fl3ll . In fact 
considering Eq. 0, in such a limit one obtains A + w 2, so 
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Eq. (0 is one of the results of the present paper. 

For n ^> ksT Eq. (0 reduces to the Lindhard function 
of the two dimensional electron gas (2DEG) while 
for /i = the prefactor before the square brackets is equal to 
-(4\n2/Tr)k B T/h 2 v F . Thus, as anticipated, for fi = k B T 
( j the only process that contributes to nonzero Imx is due to 
Eq. 0. 

Comparing Eqs. and 0, we notice that the domains 
where the function lmx±(uJ, K), which 'counts' the num- 
ber of particle-hole excitations of energy u> and momentum 
K, is nonzero are different: |w| > v F K for Eq. and 
| a; | < v F K for Eq. @. At the formal level this difference 
in domains where particle-hole excitations are allowed comes 
from the fact that Eq. describes processes where two in- 
volved quasiparticles are from the same branch of the spec- 
trum (see Eq. 0), while Eq. 0, as mentioned above, origi- 
nates from processes where two quasiparticles belong to the 
different branches of the spectrum (see Eq. 0). Using the 
Dirac equation analogies, one can say that for nonzero elec- 
tron mass A, Eq. would describe electron-electron scatter- 
ing, while Eq. is related to the process of electron-positron 
pair creation. 

As predicted in Refs. 1 14] and experimentally confirmed 
(see e.g. Ref. 0]) for T > 0.2T F (where T F is the Fermi 
temperature of the electron gas), in metallic electron gases, the 
coupled acoustic and optic plasmon modes occurring in the 
layered system do contribute to the transresistivity resulting 
in its substantial enhancement. 

In the case of graphene we do expect such enhance- 
ment, since the semimetallic band structure of graphite allows 
nonzero lmx(oJ,K) for \oj\ > v F K even for T = 0. Due 
to this peculiarity, and at variance with usual 2DEG Coulomb 
drag, the transresistivity plasmon enhancement occurs already 
at low temperatures, T > (see Fig. \l\ inset (c) where the 
transresistivity is plotted for low values of the temperature 
rescaled by the characteristic Coulomb energy, see below). 

Following the method described in Ref. |6], we have de- 
rived the plasmon dispersion relations for small K and ui 
(\u)\ > v F K). Their expressions are given by 



kw ac = 2Vln2\ — k B TKd and (8) 
V sd 



hui opt = 2yf2-Sh^2\ —k B Ty/Kd (9) 
V ed 

Notice that, apart from the factor vln 2, Eqs. and are 
formally similar to the well known electron liquid expres- 
sions, except that the Fermi energy /i (/i = in the undoped 
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system we are considering) is substituted by the characteris- 
tic Coulomb energy Ec = e 2 /(ed). Therefore the system 
behavior will be in this case characterized by the crossover 
between Coulomb and thermal energy. This allows us to de- 
fine the dimensionless temperature T a = ksT/Ec such that 
for T a > > 1 the system will be in a non-degenerate regime, 
while for T a < < 1 the system will be dominated by Coulomb 
interactions. This simple energy scale-related analysis con- 
firms that, even for small temperature (T a « 1), plasmon en- 
hancement will be important in graphene systems (see Fig.Q 
inset (c)). 

If we introduce the dimensionless quantities Co = 
hu/^E c k B T, K = Kd, we notice that Eqs. © and © ac- 
quire the universal form 



(10) 
(11) 



which is plotted in Fig.^ inset (a). If in addition we define the 
"fine-structure" constant a -1 = (hvF/d)/Ec = hvFt/e 2 , 
we find an important scaling property of the transresistivity in 
respect to the interlayer distance d: Eq. ([0 in fact acquires the 
form 
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is a universal function of T a and a, independent of the inter- 
layer distance d. 

In Fig. [2 we plot pu as a function of T a . The solid line 
corresponds to use in Eq. Jl 3i the full T-dependent response 
function x(T) (sum of Eqs. (|6j and (0); the dashed line cor- 
responds instead to the approximation Eq. (|6j. The effect of 
including x+> i- e - including temperature dependent effects, 
is quite significant. At variance with the usual Coulomb drag, 
the strong transresistivity enhancement is due in this case both 
to the enlargement of the single-particle phase-space - which 
includes now intra-band excitation - and to collective optical 
and acoustic plasmons. All these excitations are forbidden at 
T=0 and not included in \- ■ 

The inset (b) of Fig. ^presents the calculation of p2i(T) 
from Eq. Q in respect to temperature, for two different in- 
terlayer distance (d = 375A and d = 1000A) and a typical 
2DEG density. This clearly shows that the dependence on d 
of the transresistivity is very strong. In this respect it is even 
more remarkable that both curves of inset (b) collapse onto 
the solid-line curve of the main panel when the reduced units 
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FIG. 1 : Universal transresistivity pu vs. dimensionless temperature 
T a . The dashed line corresponds to the approximated calculation 
using \ = X-- Inset (a): rescaled universal plasmon dispersion 
relation Hi vs. K. Inset (b): transresistivity pn vs. temperature for 
two different value of the interlayer distance d (as labelled). Inset 
(c): as main panel, but for low values of T a In all calculations we 
have used Vf = 9.7 x 10 5 m/s ande = 6. H 



are considered. It would be very interesting to check exper- 
imentally such scaling behavior, which decouples the effect 
of layers' separation from the interlayer Coulomb interaction. 
Eq. dl 2l i in fact suggests that the main effect of separating two 
interacting graphene layers by a distance d is to renormalize 
the layer electron densities according to such a distance. We 
emphasize that, when considering usual Coulomb drag exper- 
iments between quantum wells, such effect would be hidden 
by unavoidable experimental fluctuations of different param- 
eters, such as finite quantum well thicknesses or doping. To 
check our predictions, we suggest a Coulomb drag experiment 
between two large radius carbon nanotubes, coaxial 1 1611 or 
with parallel axes. In this way the aforementioned fluctuations 
would be automatically avoided allowing the experimental ob- 
servation of such a scaling property. In the remaining part of 
the paper, we concentrate on the plasmon behavior and on the 
effect of this scaling property on plasmons. In Fig. [2] we plot 
the rescaled integrand 87r 2 /(o), K; a) j \u\fT~ a ) (see Eq. dl3> ) 
as a function of Co for different K and T a . Let us first focus on 
Fig-Etd). It shows the low K regime, in which the approxi- 
mate expressions Eqs. dlOt and Jl 1> describing the plasmons 
hold. The dimensionless temperature T a is varied by a factor 
10 between T a = 0.1 and T a = 1000, as labelled O. The 
vertical dashed-double-dot line corresponds to the position of 
the optical plasmon (from Eq. dl 1» . and the dashed line to the 
acoustic plasmon one (from Eq. dlOt . Notice that plasmons 
are present in the domain \u)\ > hvp only, so if the value of 
uj corresponding to the plasmon does not belong to such a re- 
gion, the corresponding plasmon disappears. This is the case 
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tally observable. Most of the calculations presented were 
made for the simplest gapless form of the quasiparticle spec- 
trum in graphene, as expected from tight-binding calculations 
in the absence of interactions. These results, can be easily 
generalized to take into account the possible opening of a di- 
electric gap A in the quasiparticle spectrum 1 1, 4]. The inter- 
play of this gap with the temperature and chemical potential 
may lead to even richer physics and applications. 
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FIG. 2: Transresistivity integrand 8tv 2 I(u), K) / (aVTa) vs Q for dif- 
ferent values of K (as labelled). Panels (a), (b), (c), (d) and (f): 
Five values of T a (separated by a factor 10) have been considered, 
T a = 0.1 — > 1000, each corresponding to a different solid line (see 
labelling in panel(d)). Dashed (dashed-double-dot) line corresponds 
to analytical position of the acoustic (optical) plasmon. Panel(e): 
Only two values of T a are considered (as labelled); dashed lines cor- 
respond to calculations by using the approximation \ = 

for the acoustic plasmon at T a = 0.1 and T a = 1. Fig.|3d) 
clearly shows that the scaling law Eqs. (II Oi and il Tj are well 
satisfied, irrespectively of the value of T a . In addition, due 
to the rescaling 8ir 2 / {a^/T^) , the strength of the optical plas- 
mon becomes independent of T a . 

By inspecting Fig.|2ja), (b) and (c), we see that the behav- 
ior predicted by Eqs. (II 01 and ( II II is qualitatively respected 
even for K of the order of unity, especially as long as the 
optical plasmon is considered. In particular these equations 
predict an overlap of the two plasmons at K = 2 (dashed line, 
panel(b)). Indeed the two plasmons join for a reasonably close 
value of K (see panel(c), K — 4 ). For increasing K (large K 
regime) though, a single plasmon-structure is found, i.e. the 
plasmon behavior differs even qualitatively from Eqs. (1101 and 
(O (see panel (f), K = 10). 

Fig.|2je) compares 8tt 2 I(uj, K; a) / (a\/J\) calculated us- 
ing x — X- + X+ (solid line), to the approximation obtained 
using Eq. (|6) (dashed lines), for two values of T a . As can 
be clearly seen, the integrand is influenced not only by the 
plasmon structure, which enhances it for \uj\ > hvp, but 
also by the single particle intraband excitation continuum: for 
\lu\ < hvp, the integrand would be otherwise vanishing (see 
dashed lines). 

In this Letter we have presented the first, to the best of 
our knowledge, calculation of Coulomb drag effects between 
graphene layers, which include temperature-dependent ef- 
fects. We have also predicted a universal behavior of the 
Coulomb transresistivity, which is interlayer-distance inde- 
pendent and suggested that this behavior may be experimen- 



* Electronic address: damico@isiosf.isi.it 
' Electronic address: sharapov@bitp.kiev.ua 
[1] Y. Kopelevich, P. Esquinazi, J.H.S. Torres, R.R. da Silva, 
H. Kempa, F. Mrowka, and R. Ocana, Studies of High 
Temperature Superconductors, Vol. 45, pp. 59-106 (2003); 
cond-mat/0209442 
[2] J. Gonzalez, F. Guinea, M.A.H. Vozmediano, Phys. Rev.B 59, 
2474 (1999). 

[3] J. Gonzalez, F. Guinea, M.A.H. Vozmediano, Phys. Rev.B 63, 
134421 (2001). 

[4] D.V Khveshchenko, Phys. Rev. Lett. 87, 246802 (2001); 
E.V Gorbar, V.P. Gusynin, V.A. Miransky, and LA. Shovkovy, 
Phys . Rev. B 66, 045 1 08 (2002) . 

[5] M.B. Pogrebinskii, Fiz. Tekh. Poluprovodn. 11, 637 (1997) 
[Sov. Phys. Semicond. 11, 372 (1977)]; P.J. Price, Physica (Am- 
sterdam) 117B, 750 (1983). 

[6] See e.g. A.G. Rojo, J. Phys.:Condens. Matter 11, R31 (1999), 
and references therein. 

[7] A. Kamenev and Y. Oreg, Phys. Rev. B 52, 7516 (1995); 
K. Flensberg, B. Y.-K. Hu, A.-P. Jauho, J.M. Kinaret, 
Phys. Rev. B 52, 14761 (1995). 

[8] I. D'Amico, G. Vignale, Phys. Rev.B 68, 045307 (2003). 

[9] A. Larkin and A. Varlamov, Theory of fluctuations in supercon- 
ductors, Clarendon Press, Oxford, 2004. 
[10] I.A. Luk'yanchuk and Y. Kopelevich, cond-mat/0402058 
[11] A.M. Tsvelik, Quantum Field Theory in Condensed Matter 

Physics, Cambridge University Press, 1995. Chap. 12. 
[12] R. Pisarski, Phys. Rev. D 29, 2423 (1984); T.W. Ap- 
pelquist, M. Bowick, D. Karabali, and L.C.R. Wijewardhana, 
Phys. Rev. D 33, 3704 (1984). 
[13] By inspecting Eq. Q, we then expect that our results for the 
drag transresistivity will be more reliable for relatively small 
temperatures (hui < ksT) and finite interlayer distance d (K < 

[14] K. Flensberg and B. Y.-K. Hu, Phys. Rev. Lett. 73, 3572 (1994); 
Phys. Rev. B 52, 14796 (1995). 

[15] We underline that, depending on the different d considered, the 
same dimensionless temperature may correspond to low or high 
temperatures T - for example for d = 1000A and d — 375A, 
T a = 1 corresponds to T = 28K and T = 74K respectively. 

[16] A.M. Lunde and A.-P.Jauho, Semicond. Sci. Technol. 19, 1 
(2004). 



